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CHAPTER I 


INTRODUCTION 

Three-dimensional turbulent boundary layers occur in both 
external and internal flows. Examples are boundary layers on swept 
wings and boundary layers on turbomachinery endwalls between 
compressor blades or turbine blades. Gaining physical understand- 
ing of such flows is very important for both analysis and design pur- 
poses. 

The present work treats turbomachinery endwall boundary 
layers, pointing out some basic physics of such flows and suggest- 
ing a model which attempts to describe parts of the flow-field 
analytically. 

Figure 1 shows a schematic of a turbomachinery flow-field, 
indicating the regions treated by the present work. Region I is the 
part of the endwall flow which is not influenced by the blades and 
which is pressure-driven, the crosswise pressure gradient creating 
the secondary flows which give this part of the flow its three- 
dimensional character. Region II is the blade/endwall juncture re- 
gion. where the three-dlmenslonallty of the flow is created by 
secondary flows induced by interference effects due to the mutual 
interaction of the two boundary layers on each side of the juncture. 
There is another juncture region between Region I and the suction 


side of the other blade, which Is omitted from Fig. 1 for clarity. 
The cross-hatched plane depicts a typical streamwlse station in 
which the present analysis Is done. 

Region I Is treated much more extensively In the literature 
than Region II. Coles HI was apparently the first to suggest the 
possibility of treating three-dimensional turbulent boundary layers in 
vector form, whereby the entire velocity profile is described by the 
sum of two vectors: a wall component and a wake component. 

Based on the experimental data available to him. he tentatively stat- 
ed that If the wall component Is in the direction of the surface shear 
stress, then the wake component would be very nearly parallel to 
the direction of the pressure gradient. Johnston [2] suggested the 
approach of treating the secondary (or crosswise) component of a 
three-dimensional turbulent boundary layer as a function of the main 
(or streamwlse) component and of certain additional parameters. 
He was thus led through examination of data to the proposition of 
his well known triangular model to describe the form of the velocity 
hodograph. Based on his own experiments and those of others, he 
suggested two explicit relationships to describe the triangular model: 
one for the region adjacent to the wall and the other for that near 
the freestream. His model is successful enough In the regions 
away from the vertex of the triangle but fails to match the data ade- 
quately in the vicinity of the vertex. As will be shown later, this 
region corresponds to the transitional region from the wall- to the 
wake-sublayer. in spite of its relative success, Johnston's model 


2 



suffers from being a product of data-correlatlon without the benefit 
of a supporting analysis. 

The main phenomenon In juncture flows is the corner vortex 
[3-61, which Is composed of two sources: the ’horseshoe’ vortex, 
created by the wrapping of the oncoming vortlclty around the blades, 
and the ’passage* vortex, which rolls up from the skewed boundary 
layer-induced crossflow. The blade leading edge region is predom- 
inantly Invlscld and rotational. Shabaka's [51 experiments Indicate 
the existence of a Reynolds stress-induced ’boundary jet’ outflow 
from the juncture over the endwall. It Is Important to realize that, 
since the eddy viscosity scales are different for the blade and 
endwall. a single eddy viscosity model cannot work for juncture 
flows. 

In order to gain further Insight Into the physics of three- 
dimensional, endwall turbulent boundary layers, it became clear that 
more analysis would be helpful. To this end. Mellor's [71 two- 

dimensional asymptotic analysis for large Reynolds number turbulent 
boundary layers has been extended to three dimensions In the 
present work. The results of this analysis provide the scaling of 
three-dimensional boundary layers of both juncture - and pressure- 
driven - types, including rotational effects. The analysis of the 
former type Includes both symmetric and asymmetric cases. The 
analysis of the latter type leads to the proposition of a composite 
three-dimensional Law-of— the- Wall and Law-of-the-Wake model, util- 
izing a vector approach In the wake sub-layer. This model is com- 


3 



pared with van den Berg and Elsenaar's 18] infinite swept-wing flow 
data and Mueller's 110] forced turning flow data with good agree- 
ment. It should be mentioned that Johnston's triangle emerges as a 
particular limit of this model. 

The analysis presented here does not Include solutions of the 
equations of motion found for the various sub-regions, and therefore 
It comes short of predicting forces and moments acting on the 
blades of a turbomachinery component. 
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CHAPTER II 


ASYMPTOTIC ANALYSIS 


11 . 1 . The Pressure-Driven Case (Region I) 

A typical skewed boundary layer and basic nomenclature are 
shown In Fig. 2. For incompressible, steady flow, with no curva- 
tures, the mean part of a turbulent flow is described by the follow- 
ing equations: 

Mass conservation: 

au f /ax.= o <u 

Momentum: 

u. [au/aXj] = - ap/ax f + 

a/ax / (e 2 r // t€ 2 eau / /ax / ) <2) 

where 

u rVV ' >=pV <’ u 1- V - "/’"V"* 2 - 


5 


x =x /i , e=u*/u a . e 2 e=v/( J /)=fl _1 
II e e e f 

In the above equations, e 2 is the scale for the Reynolds stresses. 
* 

with u being an as of yet undetermined perturbation velocity scale. 
e 2 € Is the scale for the viscous stresses, which Is also given by 
. Thus e Is the relative magnitude of viscous- to Reynolds- 
stresses. 

In addition, the following boundary conditions are required: 
tf / (X,0.Z)=0 

and. assuming a uniform external flow. 

0 i x - z ) = ft} 

(for a fixed (X.Z) location ) 

Next, an asymptotic analysis is carried out (following Mellor 
[7]) for the above equations, using the following expansions: 

A. Outer Expansion 

0=Q-j +eQg+e^Qg + - • • 

7^=0 ( This assumes no free -stream turbulence ) 
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7 


and pressure: 


n expression of the form 
at a given (X,Z) location. 



U y =ee 


\ t€ u f 

y l y 2 y 3 


A A 2 A 

U=u + eu + c u + 

1 2 Z 3 


P 


A 



A 2 A 

+ ep 2 + e p 3 + • • • 


r // = 


+ ef„ 


2? 
e f 0 


'// 


’// 


J // 


These expansions are substituted into Eqs. (1) and (2) and the 
analysis proceeds in the manner of matched asymptotic expansions, 
tho details of which aro described in Appendices A. B and C. It 


should be noted that this analysis assumes that both e and e are 
much smaller than unity. The features of the resulting equations in 
the throe regions are summarized in Table 1 to the first three ord- 
ers. Tho main features of the flow, as indicated in the table, are 
the following: the outer layer is inviscid; the middle layer is influ- 
enced by Reynolds stresses to orders 2 and 3. where the 2nd order 
equations are of the boundary layer typo while the 3rd order ones 
aro no avoragod Navior Stokes equations; the inner layer is dom- 
inated by constant total shear (viscous- plus Reynolds-stresses) 
across it. with no pressure gradient or inertia effects through order 
3. The important information for the subsequent analysis is con- 
tained in the 1st ordor outer region equations and in the 2nd order 
middle- and inner- region equations. 
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As a result of the asymptotic analysis and the matching pro- 
cedure, the following conclusions emerge: 

1. The flow possesses a 3-layer structure; 

2. The Inner and middle layers have a common loga- 
rithmic overlap; 

3. The Inner layer is exponentially thin compared to the 


middle layer: e = O 



as e -• 0: 


4. The pressure gradient does not influence the Inner 
layer, hence the wall shear stress is the only driving force 
there. It is reasonable to assume that the shear field main- 
tains its two-dimensionality within this sub-layer. Thus the 
flow is two-dimensional In the inner layer and the velocity 
vectors are co-planar in the r -direction. 

5. Since the outer flow Is In the freestream direction, all 

the turning from the r - direction to the 0 - direction 

it e 

must take place In the middle layer. 

While the first three conclusions apply also to two-dimensional 
boundary layerss, the last two are special to the three-dimensional 
case. 

See Appendix G for the effects of rotation on the equations of 
motion In the pressure-driven case. 


11.2 . The Case of Juncture Flow (Region II) 
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Figure 3 describes the main features of the flow field In a 
juncture. An oncoming endwall boundary layer interacts with a 
blade situated normal to the endwall. This interaction gives rise to 
several phenomena: the vorticlty carried by the oncoming b. I. wraps 
around the blade, thereby creating a "horseshoe" vortex; the stag- 
nation line at the blade leading edge creates a pressure gradient 
that induces skewing of the oncoming b.l. and the formation of an 
endwall crossflow, which later rolls up into a so-called "passage 
vortex"; this and the "horseshoe* vortex augment each other and 
create the corner vortex which stretches downstream without dimin- 
ishing its strength for quite a distance downstream of the blade 
leading edge; If the latter is blunt, the oncoming flow may separate 
upstream of it. creating one of several possible three-dimensional 
separation patterns, such as a saddle-point or a node. 

As the oncoming flow hits the blade, a boundary layer 
develops over it. This b. I. interacts with the one which has been 
developing over the endwall. resulting In a corner (or juncture) flow 
which is investigated in this section. 

11-2.1. Scaling and normalized equations 

The flow variables in the juncture are normalized as follows: 

u l - VV p ■ P/ P °l- T„ = 
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X = x'/x ref . y = y/\ w a. z = Z/\ b 6 (X = ^27c~ f ) 

where (x.y.Z) is a locally cartesian coordinate system. 

The scaling of Y and Z Is based on the results of the 
analysis of the three-dimensional boundary layer outside the juncture 
region (Region I). which Is performed In Section 2.1. (See also 
Appendix D) . When the equations of motion are written In terms of 
these variables, the factors X^/fX^Alsi and 6/(X^A) = C ap- 
pear in them (see Appendix FI) and these require additional con- 
sideration. By using the (l/n)th power law for turbulent flat plate 
flow as a model, it can be shown (Appendix FI) that the equations 
of motion for the juncture flow become 


du av i aw m o 
ax T e ay ( a z 




a p . a r.-i ,2i au 
e w e w ax 


= _ + 

ax ax I* 


t e 


2 7 1 
W XX 1 


a r 2 ^ bu 

2 ^ 1 


L a 

r. 

-12 a ay 

+ € ayLVw ay 

* T yx\ 

+ 

C az 

i c 

e w e w az 

av > w a\/ * 

u/ av 


ap A 

a 

r,-i 2 a 

ax * v ay c 


£ 

ay + 

ax 

L* e w € w 


+ e 


2 7 1 

w xy J 


(3) 
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4. t ±r 2 s dv t 2 r i 3 2a av 2 -i 

* ayl e w € vy ay +e w r yyJ + C az K € w € w 3Z + e w T zy\ 

•-*!?♦£ [«-' «!&,« 

♦ < 2o 4 »■„] * 

+ * J_ r e 2 - avv .20 2 T 1 iAL-1 2j aw 

5 ay L € w e w ay + c e w T yzJ + £ az K e w e w az 

+ < 2o e 2 t 1 
w zzl 

where, in analogy with the treatment of Region I. 

V» “ •"“Vw"- W '°e 

w 

and the additional notation is 

C s ( W 1-2 °‘ 0 s 1/ <l-2").€ ■ x r0f /(X w A) 
The boundary conditions are: 


O = V = IV 


= 09 fy = o.z>o 
[z=o.y>o 


11.2.2. Asymptotic treatment 

Two Independent parameters appear in Eqn. (3). namely e 

w 

and C (the 3rd parameter, e is related to e as shown in Ap- 
pendix B) . While Is a small parameter everywhere along the 
juncture. C increases from 0 at the blade leading edge to 0(1) at 
the trailing edge, and may thus be viewed as a measure of both 
streamwise location and extent of asymmetry. For the purpose of 
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asymptotic analysis of Eqn. (3), a new parameter, €^, is used as 

A 

a basic perturbation parameter, to be determined (along with e^) 
from matching conditions (Appendix F2) . in order to avoid double 
expansions in terms of € w and C* It is intended to stretch the junc- 
ture coordinates Y and Z in such a manner that the juncture 
c c 

equations derived from Eqn. (3) are symmetric w. r.t. these coordi- 
nates as well as w. r.t. V and VV regardless of whether the flow 

c c 

itself is symmetric or not. Such a coordinate transformation leaves 

e as the only perturbation parameter, enabling single expansions as 
c 

for the Region i treatment. The details of the coordinate stretchings 
are given in Appendix FI, resulting In the following expansions, valid 
for an arbitrary (non-zero) C’ 

Middle region 

y = y/<*€ c ). z = c z/((c c ). € = 4 v (€ c®c > 

u “ u i + e c u 2 + e c U 3 + * “ 
v - € C h + ‘ e c v 2 + e c v 3 + * * * ) 

W = e c [ W ^ + e c"2 + e c w 3 + * * * ) 

p " P] + e c p 2 + € c P 3 + *•* 
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T II • (W)* [' + («' 2 ° - ' ) • 


'// c 2 // c 3 // 


fO, / 1 4 ; 

*°/* = ll./ = 


*z 

z 


Inner region 


y = y/(«e c e c ). z = c Z/((e c € c ) 


t . A a p A 

U = "l ’ Vs ♦ e c°3 * 


v = Vc ( v "i + Vs * V 3 + • ) 

w = Vc (»1 + Vs + V 3 + ■ • ■ j 


P " P 1 + Vs + e c P 3 * ' ' ' 


r // 55 [«</*»]* [' + (i" 2 ° - i)o /2 ] 


[ 


. A 2 * 

+ e n? + Oa + 

W c2 // c3 // 




As shown in Appendix F2. 


" m € c = V ,lme c = V lim « e = V ,,m K “ *f> 

* ; _ 

c 7 c 

These expansions are inserted into Eqn. (3) and the resulting equa- 
tions are sorted out to orders 1 . e c and e* for the middle and 
inner regions, as shown in Appendix FI. it is important to realize 
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that the Reynolds stresses scale differently In the endwall- and 
blade- boundary layers (Regions I). Thus, In the endwall b. I. . 


T w = f l + e w f 2 + 
W ll 'll II 


with 


related expansion for the Inner layer. In the blade b.l. : 


T. = f, + e. + 

b il V b2 n 


and similarly for the Inner layer. The Juncture region contains Rey- 
nolds stress contributions from both boundary layers adjacent to It, 
and this is taken into account by the series expansion representation 


of T as demonstrated below: 

c ll 


llm T 

C,i 

z c — II 


1 + e w f 2 + 

'll W II 


. 1 * 2 . 


F’ 


C 2 °l'l + e w f 2 + 

'll W II 


. 1=2 


llm T 


'll 


‘‘°[\ ’ A/ *•"]■' 


* z 


U + e b f 2 + 

'll D * II 


. I = 2 


Thus, the factors £ ±So compensate for the two different scales for 
T on endwall and blade, enabling smooth matching of the Juncture 
region Reynolds stresses with those on the two Region I boundary 

layers. 
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The equations of -motion are derived In Appendix FI. while 
Appendix F2 includes a detailed description of the procedure to 
match the corner region with the boundary layers outside the 
corner. Since the equations are symmetric, it suffices to perform 
the matching to half of the corner, say the sector between the 2- 
axls and the bisector of the right angle between this and the Y-axis. 
The resulting matching lines are indicated in Fig. 4 by the numbers 

1. 2 and 3. Important results from the matching procedure are 
given In Appendix F2 and summarized below: 

1. The pressure in the corner region Is Imposed by the 
external flow; the gradients of the external pressure normal 
to the solid boundaries vanish along the corner. 

2. The streamwlse velocity has a common logarithmic 
overlap between the corner middle region and the inner 
layers on both endwall and blade. 

3. The streamwlse velocity component In the inner region 
vanishes to 0(1). 

The final version of the equations of motion is given in Table 

2. These equations indicate the following main features of turbulent 
corner flows: 

1. The pressure field In the corner region Is Imposed by 
the outer flow. 

2. In the middle region, the flow is inviscid and without 
Reynolds stresses to 0(1). and is described by an Euler- 
type equation for the streamwlse momentum. This same 
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momentum equation becomes of a parabolized Navier-Stokes 

2 

type with Reynolds stresses only, to 0(e). The OCe ) flow 
is described by a Re-averaged N-S-type equation for the 
streamwlse momentum and by parabolized N-S equations for 
the two crosswise momenta. Again, these equations in- 
cludo Reynolds stresses only. 

3. The inner region is characterized by the fact that the 
crosswise derivatives of the total shear (viscous + Reynolds 
stresses) on each wall are equal In magnitude and opposite 

2 

in sign. Also, the 0(e ) flow Includes pressure gradients 
normal to both walls, which are balanced by crosswise 
Reynolds stress gradients. There are no inertia effects in 

2 

the inner region through 0(e ) at least. 

It is noted that, unlike the Region I equations, those of Re- 
gion II cannot be solved starting at the lowest order and continuing 
to higher orders based on previous solutions. Rather, one has to 
use symultaneously equations belonging to different orders, to have 
enough non-trivial equations for the number of unknowns [15,161. 
Specifically, in the middle region the continuity and streamwlse 

momentum equations of 0(e m ) have to be coupled with the two 

crossflow momentum equations of 0(e™ +2 ) to solve for the 0(e m ) 

velocity components and 0(e™ + ^) pressure; In the inner region the 

continuity and streamwlse mometum equations of OCe™ > must be 
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solved symultaneously with the crosswise momentum equations of 
/ 7 ? i*2 A 2 

O ( e e ) to get a complete flow picture. This means that con- 
siderably more terms of the asymptotic series must be retained In 
Region II than In Region I. especially In the Inner region. Note 
also that a 'composite profile* does not exist for the juncture flow, 
since the outer, middle and Inner regions do not share common 
surfaces and thus do not have common asymptotes (see Fig. 4) . 

Finally, the extent of the juncture region for an arbitrary 0/A 
ratio is of interest. This is determined by letting the coordinates of 
the various subregions be of the order of their scales. Thus 


y c = 0(1) - Y = O 


(fe c ) - Y = ofej e c £A) = 0 [*Au / u 

c w 


K = on> - r = o( £ e o aj- y = o(e; 


1 A , 
€ c e c< 


1 - 0 

Av 

J - 0 

* A u 
c r 

w 


Z c =0(1) - Z 


- ° 

rU 


T b i 


z, = 0(1) 


2=0 ( e c a c €/c ) - 2 = 0 [ € ~b € ch a ^ /c ] 


= 0 


L JL v 

i : a , u 

c T b 


Figure 4 shows the various sub-regions for a juncture flow with 
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nonzero but otherwise arbitrary 0/A ratio. The corner vortex, dis- 
cussed earlier in this Chapter, is situated in the Juncture middle re- 
gion. 

For rotational effects on the equations of motion In the junc- 
ture, consult Appendix G. 
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CHAPTER III 


PROFILE RELATIONS FOR PRESSURE-DRIVEN BOUNDARY LAYERS 

111-1- General Discussion 

Based on the results of the pressure-driven 3D analysis, the 
general forms for the Law-of-the-Wall and Law-of-the-Wake are the 
following: 

Law-of-the-Wall 



Law-of-the-Wake 



where 
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V 


Y 
6 

Matching these two laws yields the logarithmic overlap (see Appendix 
D) : 

f(x.t).z) " -K _1 Int) + A '.u-0 

u + (y + ) " k 1 /ny + +B.y + 

where 

K2 .41,8 a 5.0. A' = 2k 1 II 
which are Coles' parameters. 

If u* Is now chosen as u then the relationships between F 

T 

and i/ + on the one hand and the 2nd order middle- and Inner- 
variables on the other, become (see Appendix D) : 


1/2 



This enables the construction of the composite velocity profile, fol- 
lowing the general scheme 

fcompos/fe] _ (outer profile) -(common outer-middle 
[ profile J 

asymptote) +( middle profile)- (common middle-inner 

asymptote) + ( Inner profile) 
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Written out explicitly, the composite profile, to order e. has the 
form (see Appendices C and D) : 



and y is the angle between F and 0 . In order to be able to 

© 

use this profile in actual calculations. It Is necessary to provide ex- 
plicit models for the Wall- and Wake-Laws. 

III. 2. LAW-OF-THE WALL 
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Based on the results of the asymptotic analysis, the 
Wall-Law Is written In the form: 


III. 3. 


= <(/).* - < 4 ' 

with components 

u* = u + cos$ w . u + = u + slntf^. 

For f Musker's C121 explicit expression Is chosen: 

f </ + > = 5.424 tan 1 j 2y -,0 7 " 1 5 } 


+ log 1Q 


(y + + 10 . 6 ) 9 ‘ 6 

+ 2 + 
y -8. 15y + 86 


-3.52 


LAW-OF-TH E-WAKE 

In vie * of the conclusions from the asymptotic analysis, 
the wake function Is defined In the form: 



T 


with streamwlse component 
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F x = Fcosy = 
and crossflow 



u 


T 

component 


U 


F = - Fslny 

z u 


See Fig. 5 for basic nomenclature. 

The two-dimensional defect function may be expressed by 
Musker's (12] formula 


F 2D = * 1 (2Il[l-37? 2 + 2t? 3 ) - ln V - 7> 2 (l-77]] 


(5) 


where the Wake Parameter, n. Is determinable from 1 (see Appendix 
El) . 


211 - In 


i+rr 


k Re 


+ K 


( S_ M - 


= 0 


with 


°« ' 1 C' 0 *J ay ' 

and 



X 

(oos^) ,/2 


There Is not an unique generalization of F 2£> to three dimensions. 
Possibilities for Fg^ Include: 
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1 . 


u e - u 


/ u 





u . 

T 


The first possibility was rejected since it led to poor data - correla- 
tion; the second is difficult to test or use; the third possibilty is the 
simplest and is the one used henceforth. Thus 


F 


' 2D 
cosy 


(6) 


The geometry of Fig. 5 provides the Wake-Law formulation: 


u = 


1-2X V cosy+ 


tx-V] 


1/2 


0=tan 


-l r x V 

li- X- 1 


^sinyl 
F cosr * 


(7) 


or, equivalently, 

ta,vy = { F coiT" _1 } tan0 <8> 

In order to complete the formulation of the Law-of-the-Wake, 
It Is necessary to determine the function y<7?) . Figure 6 shows a 

plot of experimental values of r s where y m is the maximum 

y m 

value of y. The data are from Refs. 8 and 10. It Is observed that 
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the lower portion of the plot is different for different stations and 
that these Individual branches merge Into a single curve further up. 
This suggests that r could be describable by a combination of two 
parts, one of which Is a function of 17 and of some physical param- 
eters peculiar to each measurement station, while the other Is a 
function of 17 only. Pursuing this direction, one arrives at the for- 
mulation described below, the details of which are worked out in 
Appendix E2. In the logarithmic overlap region, y is determined 
from the solution of 


2 

dy _ Xcos y 

drj ~ p 2 

*2 D tany 


U 

K-n 


2 D 


Xcos^y 


-1 



( 9 ) 


where U Is given by Eq. ( 7 ) and F 2£) by Eq. ( 5 ). 

The initial condition is 

y - r 0 « V = T ) 0 (9a) 


where 77^ is the point where the Wake— Law crosses 0^ (see Appen- 
dix E 2 . and Figs. 7 and 8) . 

In the defect layer, y Is chosen In the form 

r - ( 1 -77 ] “ ( 10 ) 

m 


as shown In Appendix E 2 . where 
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(10a) 


a = 



fill 
ld7j Jt? 1 



. 0 = 


(’-M 


with subscripts 1 denoting the point of switching from the portion of 
the curve described by Eq. (9) to the one described by Eq. (10). 
(see Appendix E2) . The dimensionless distance 7?., is a function of 
Reynolds number and pressure gradient. A typical value Is 0.035. 
Figure 8 depicts the relation of 77 1 to v 0 - and Appendix E3 outlines 
the procedure to determine T)^. 

At the upper edge of the defect layer y attains Its maximum 

2 

value, y . which Is calculable from 
m 


m 


tan 


-1 




e 


HD 


where 0 Is the local freestream direction relative to the direction at 
Q 

the upstream position. Eqn. (11) is a result of integrating a 
shearless version of the boundary-layer equations along main-flow 
streamlines. This approach Is based on the assumption that a 
shearless analysis can yield valid results In the outer part of 
boundary-layers which are well developed before being perturbed by 
a main-flow turning. Under such circumstances crossflows develop 
mainly as a result of the reorientation of vorticlty In the layer, while 
frictional effects are felt only near the wall. The determination of 
y thus requires knowledge of the external flow field only. Figure 9 
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depicts a typical comparison between the above model for y and ex- 
perimental data. 

It should be noted that the entire analysis described above did 
not Involve the actual solution of the equations of motion. There- 
fore, four parameters are required as input In order to enable cal- 
dilations. These are: 

1. c, 

2 - K 

S. Re 6 

4 *e 

In addition, the determination of y m requires knowledge of the 
external flow field. In a real situation these would be unknowns, 
and the equations of motion would have to be solved, subject to a 
prescribed pressure distribution and initial velocity profile. However, 
regardless of the method of solution, the above model could then be 
incorporated Into the solution scheme, since it does not Introduce 
any additional unknowns of its own. 


Ill -4. Comparison with Experimental Data 

The data of Refs. 8 and 10 have been used to check the 
validity of the proposed models for the Wall- and Wake-Laws. Both 
sources of data use flow geometries which lead to a build-up of ad- 
verse pressure gradient with eventual three-dimensional separation. 
The measurement stations chosen for comparison between theory and 
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data were all fully attached. Results of two stations from Ref. 8 
and three from Ref. 10 are reported here: In station C4 of Ref. 10 
an 8% crossflow Is present: In station D5 the crossflow Is about 

12% and In station F7 It reaches about 25% of the streamwlse flow 
component. In station 7 of Ref. 8 the crossflow Is about 22% and 
in station 8 it reaches 27%. Figures 10A through 10E show com- 
parisons between the present model and the data of the representa- 
tive measurement stations. The figures Include comparisons for both 
the hodograph and the resultant velocity profile (In terms of magni- 
tude and direction) . The agreement Is good for all stations. In 

particular, it is noted that the model is able to predict the smooth 
transition from the wall- to the wake-region. The discrepancy 
between the model and the data of Ref. 10 Involving the magnitude 
of the velocity (most pronounced in station F7) may be due to the 
fact that the experiment of Mueller has a curvature effect, which Is 
not accounted for In the model. 


III. 5. Comparison with Johnston's Triangular Model 

Johnston [2] assumed the following model for the 
hodograph: 


U 

z 

In 


contrast, the present model assumes no geometric relationship. 
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Rather, the hodograph results on the basis of the fundamental as- 
sumption that the wake function Is a vector quantity. In comparing 
the two models. It Is observed that the first leg of Johnston's trian- 
gle corresponds to the Wall-Law of the present model, while the 
second leg corresponds to a special limit of the present model's 

Wake-Law. namely, r = 1 with A = tan-y . However, this Is an 

m 

Inadequate choice since It falls to correctly describe either the tran- 
sition from the wall- to the wake-region or the magnitude and loca- 
tion of the maximum crossflow (see Fig. 7). These shortcomings 
are apparent from Johnston's own polar plots (see Ref. 2) . The 
other basic difference between the two models stems from the fact 
that whereas Johnston regards the vertex of his triangle as the edge 
of the wall sublayer, the present model locates this point at t ) q . 
corresponding to Vq (see Fig. 7). 
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CHAPTER IV 


CONCLUSIONS 

An asymptotic analysis has been carried out for three- 
dimensional, endwall turbulent boundary layers. This analysis Indi- 
cates that, as In the two-dimensional case, the flow Is multl- 
structured, possessing an Inner (or wall) layer; a middle (or de- 
fect) layer; and an outer layer, and the pressure in the b. I. Is Im- 
posed by the outer flow. 

In the blade/endwall juncture, the defect region is describable 

2 

by parabolized Navier-Stokes equations through 0(c) except for the 

o 

streamwlse momentum of 0(c) which is a Re-averaged N-S-type 
equation. Rotational effects appear In the form of centrifugal terms 

o 

to 0(1). Coriolis terms to 0(e ) and there Is no influence on the 
0(c) momentum components. Viscous stresses are absent from the 
defect region. The inner region Is characterized by the crosswise 
total shear gradients on each wall being equal in magnitude and op- 
posite in sign, and no inertia effects, except for centrifugal ac- 
celerations in the case of a rotor. 

It Is also found that the streamwise velocity has a common 
logarithmic overlap between the corner middle region and the inner 
layers on both endwall and blade. 
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The same equations which describe the non-symmetrlcal Junc- 
ture flow also describe the symmetrical case, but the coordinates 
and Reynolds stresses scale differently In the two cases. 

In the pressure-driven part of the endwall boundary layer out- 
side the juncture, the Inner layer Is dominated by constant total 
shear across It. Inertia and pressure gradients having no significant 
influence there. The middle layer is driven by both pressure gra- 
dients and Reynolds stresses, and may be described by boundary 
layer-type equations to 0(e) and Re-averaged N-S equations to 
2 

0(e). These two sub-layers share a common logarithmic overlap 
region. The outer layer Is inviscld and generally rotational, 

descrlbable by the Euler equations. 

The results of the analysis indicate that the wall sub-layer Is 
two-dimensional, flowing in the local wall shear direction, and all 
the turning from this to the freestream direction takes place in the 
defect sub-layer. 

Based on the results of the analysis for the pressure-driven 
part of the b. I. . a model is proposed for the three-dimensional 
counterpart of the Law-of-the-Wall and Law-of-the-Wake. which 
agrees quite well with experimental data. 


32 



BLADE SURFACE 
BOUNDARY 
LAYER ^ 


UNSTEADY 

UNIFORM 

VELOCITY-. 



HUB 

CORNER 

VORTEX 


END WALL 

BOUNDARY 

LAYER 


ROTATION 


FIG. 1 SCHEMATIC OF FLOW-FIELD, SHOWING REGIONS 
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FIG. 2 SKEWED BOUNDARY LAYER AND BASIC TERMINOLOGY 
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FIG. 4 ASYMPTOTIC SUB-REGIONS IN A TURBULENT 
JUNCTURE FLOW 
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U 2 



FIG. 7 HODOGRAPH WITH JOHNSTON’S TRIANGLE SUPERIMPOSED. 
ALSO INDICATES THE GEOMETRICAL SIGNIFICANCE OF 


39 



STATION: v. d. Berg No. 8 

* data MODEL 



FIG. 9 COMPARISON BETWEEN THE MODEL FOR T (^) AND EXPERIMENTAL DATA 
OF REF.0 , STATION No. 8. ( THE FIVE LOWEST POINTS BELONG 

TO THE WALL SUB-LAYER) 
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STATION: Mueller No. C4 


* + o data 

MODEL 



U. PHI/PHI W 



A 

FIG. 10 COMPARISON BETWEEN MODEL AND EXPERIMENTAL DATA OF SELECTED 

MEASUREMENT STATIONS FROM REFS. 8“ ANDiO. (SHOWN ARE RESULTANT 
VELOCITY IN MAGNITUDE AND DIRECTION, AND THE VELOCITY HODOGRAPH) 
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FIG. 11 DETERMINATION OF 
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Table 1 - Momentum equation In the three sub-layers, to orders 1. 
c. and e 2 . 



OUTER 

REGION 

MIDDLE 

REGION 

INNER 

REGION 

0 

R 

D 

J 

Euler equa- 
tions 

( No Reynolds 
stresses) 

P l = 

P 1 (X.O.Z) 
Euler equa- 
tions 

P 1 =P 1 (X.O.Z) 

A A A 

u -u =u 

4' y ' z ’ 

0 

R 

D 

€ 

U =U =U =0 

*2 y 2 z 2 
p 2 =o 2 2 

p 2 =o 

Boundary- 
Layer equa- 
tions. Rey- 

nolds 

stresses only 

P V° 

Constant total 
shear across 
sub-layer 

o 

R 

0 

e 2 

Euler equa- 
tions 

Re-averaged 
Navler-Stokes 
equations for 
X- and Z- 
momenta 
Reynolds 
stresses only 

A 

P 3 = 

Pg(X.0,Z) + 

A 

C^Kstant total 
shear across 
sub-layer 
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Table 2: Juncture equations for an arbitrary (non-zero) 

6/A ratio. (See text for definitions of coordinates). 
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Appendix A: Preliminary form of the equations of motion in the 

2 

various layers to orders 1 . e and e 
After substituting the asymptotic expansions for the three layers 
into Eqns. (1) and (2) and collecting terms of the same order, the 
following equations are obtained: 

1 . Outer layer 
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au 1 av 1 aiv 1 
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2. Middle layer 
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Appendix B: Matching procedure 


1 . Formal matching of outer and middle layers 
U(x.Y.Z) Is expanded by a Taylor series around Y = 0: 


dU 


1 


2 

2 a (A 


U(x.Y.Z) = U (.x. 0.Z) = e y (x.O.Z) + e 2 ;r <x. 


ay 


ay 2 


2 ^^2 

+ . . . + e U 2 (x.O.Z) + e y -gy- (x.O.Z) + . . . + 

e 2 U 3 (x.O.Z) + ... = 

2 

= u 1 (x.y.Z) + e u 2 (x.y.Z) + € t/g (x.y.Z) + . . . ; y -» 
Thus 


u i (x.y.Z) 'a, (x.O.Z): y - • 


au 1 

u 2 (x.y.Z) "y (x.O.Z) + L/g (x.O.Z) : y -* <*> 


. £ 


at/. 


t/g (X.y.Z) " -g 1 (X.O.Z) t y -gy~ (x.O.Z) t 


u 3 (X.O.Z): y - 


and similar relationships hold for the crossflow component, w. 
Following the same procedure as above, yields: 


O.Z) 


(B1) 
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\/ ] Cx .0.2) = 0 


(B2) 


VgCx.O.Z) = llm 

y - ‘ 


3V 1 

v 1 (x.y.Z) - y -gy (x.O.Z) 


( B3) 


V 3 (x.O.Z) = lim 

y - < 


av 2 

v 2 (x,y.Z) -y -gy (x.O.Z) 




P^ (x.Z) - P 1 (x.O.Z) (since dp^/dy = 0) (B4) 


ap 


P 2 (x.Z) - y -gy- (x.O.Z) + P 2 (x.O.Z); y - ® (B5) 


y 2 dP 2 

P 3 (x.y.Z) ' Y — y (x.O.Z) + y -g^ (x.O.Z) + P 3 (x.O.Z); y - » 


f (x.y.Z) "0; y - «*>, for n = 1 , 2. 3. • • • 

// 

(since no freestream turbulence was assumed to exist) 

2. Formal matching of middle and inner layers 

Again, variables are expanded in Taylor series around y = 0: 

« a 

V(x.y.Z) = e (x.O.Z) + e e y -gy (x.O.Z) + . . . + 
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e 2 v 2 (x.0.2) + 

2 A A dv 2 3 . , 

+ e e y (x. 0,2) + . . . + e v 3 (x.C 

e y -rr^ (x.0.2) + ... = 

8y 

A A A 

= e e v 1 (x. y. 2) + 

2 A A A 3 A A A 

e e Vg (x, y. 2) + e e Vg (x. y. z) + 

A 

y -* oo 

Hence 

v } (x.0.2) = v 2 (x.0.2) = Vg (x.0,2) = = 0 
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v f (x.y.2) ' y (x.0.2): y - / = 1.2.3. • ♦ ♦ 

2 2 a 9p 3 
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-» oo 


P 3 (x.y.Z) " P 3 (x.O.Z) ;y 


and slmilarfy 

t (x.y.Z) - t (x.O.Z): y - ~. n = 1.2.3. -•• 

// n IJ 

The U and W components require a special treatment. 

0i/, 

U (x.y.Z) = u 1 (x.O.Z) + e y (x.O.Z) + ••• + 

e u 2 (x.O.Z) + • ♦ • + 

+ e 2 u 3 (x.O.Z) + = u 1 (x.y.Z) + e u 2 (x.y.Z) + 

e 2 Ug (x.y.Z) + y - <*> (B8) 

From Eqn. (A7) and the boundary condition u 1 (x.O.Z) = 0, 
it follows that u 1 must be of the form 

u i (x.y.Z) = f (x ,Z )y 

which is unbounded as y - ®. Thus the first order terms in Eqn. 
( B8) cannot be matched directly. Therefore Eqn. (B8) is differen- 
tiated w. r.t. y. yielding: 
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ay 3y 


y -• aa 


3u, 3 Up „ au„ 

= -r- 1 - (x.O.Z) + e -^r (x.O.Z) + e (x.O.Z) + • • • <B9) 


Since the leading term in the Inner expansion Is of 0 



and that In the middle expansion is of 0(1). the only possibility for 
a match is that both vanish Identically. A similar argument holds 
for the W-component. Thus 


= w. 


■ o. 


and from the continuity Eq. and B.C. = 0 0 y = 0 It also fol- 

lows that v. e o. Also 


_a_ 

ay 



(x.O.Z) 


0 


and this combined with the B.C. 
u 1 and w 1 at y = 0: 




= 0 uniquely determines 

= 0 


at a given (x.Z) location. 


Since y = y/e. the only possibility for a match In the remaining 
part of Eqn. (B9) Is 
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du 2 /dy “ bg/y. au 3 /ay “ bg/y.***; y -* 0 
at/g/ay " bg/y, 3t/g/ay " bg/y. •••; y -» » 


which, after integrating gives 

Og “ bg Iny + Cg, Ug “ bg Iny + Cg , • - • ; y 0 (BIO) 


« 2 ' bg Iny + Cg.Wg " bg Iny + Cg . • • • ; y - <*> (BID 


where the b's and C's are functions of x and Z. Similar expres- 
sions hold for the w and w components. This establishes the ex- 
istence of a logarithmic overlap between the Inner and middle 
layers. Inserting these results Into (B8) yields 

(cbg + e 2 bg + ...]iny = u 1 (x.O.Z) + (ebg + e 2 bg + •••) Iny + 


* € ( C 2 - + <?[ C 3 - a s) * ••• 

or 




2 

+ e br 


+ € 


- C, 



(x.O.Z) + e (Cg - 


whence It follows that 


-cine = 0(1) 


e = 0 (e ) as e 


0 


This confirms the assumption that e = 0( e n ) which allowed a single 
parameter asymptotic expansion approach. 
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3. Results from the matching procedure 
A. From the continuity Eqn. of (A4): 




0 lEq. (S 6)1 


and using Eq. (61) . 

lim v 1 = - J 

y -♦ cx> 


lim 


0t/ 

Cdx 




-i 


au } aw 1 
ax + ~az 


cfy = 


(x.O.Z) 


= - y 


3L/ 1 aw 1 
ax + az 


( B12) 


(x.O.Z) 


Also, from the cont. eqn. of (Al) 


a^i 

ay 


(x.O.Z) = - 


au n aw 1 
ax + az 


(x.O.Z) 


Substituting this and (612) into Eqn. (63) yields 


V 2 (x.O.Z) = 0 


This, together with the fact that the 1st order outer solution satisfies 
ail the B.C.'s at Y -* <*>. so that the 2nd order equations must 
satisfy homogeneous B.C.'s there as well as at Y = 0, dictate that 
the solution to Eqn. (A2) be the trivial one: 
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0 
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Appendix C: Final form of the equations of motion and their 

matching- and boundary- conditions 


1 . Outer layer 
0(1) 



0(e) 


u = V = LV = p = 0 
2 2 2 2 U 


2 

0(e) 


au„ av~ aw„ 

+ -^r + = o 

ax ay az 


K V ! ff* slf'V V W 3) 


I 

I + V 3 I ’ *3 £]{<V V W l) 
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_ _ [a 

a 

J_] 

To 1 

[ax' 

ay 

az] 

raj 


2. Middle layer 
OH) 


3 

t^l + 

3 w p 1 



dx 

3y 

az 0 




3a^ 

3u^ 

3u^ 

3 Pl 

U 1 

■ + V, 

3x 1 

ay 1 

3 Z 

3x 


0 =- 


3^ 

ay 


0 ( 6 ) 


p 

0(e*) 


a *i 

- X 1 / 

3 "i 

3tv.j 

3p-j 

3x + V 1 

3y + *1 

az' = 

“ ~az 


du 2 

j. 


dw 2 

A 






ax 

T 

+ 

3y 

az 

- u 






a 


“a) 


3Wp 


3t/, 


3w p 

3y^ 

ax 

( U l 

+ v i 

“37 

+ V 2 

3y 

+ *1 

az + 

*2 3Z " 

u i 

BWr 

! + U 2 

3w 1 


dw 2 


aw^ 

a 


dx 

3x 

+ V 1 

3y 

+ V 2 

3y 

+ az 

(■Vs) ’ 


3(, 


jy. 


3f, 


11 


3y 


au, av_ a«r 

~Bx + "57 + ~3 F = 0 
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a 

3u a 

du 2 

au 1 

au 3 

+ v i 

3x 1 

CVJ 

i-- 

ay *3 

- r— ■ + W. 

by 1 

az 


du 


3u, 


3Pq af ‘> 


"2 az + w 3 az 


av. 


XX 


ax 


°i “57 + v 


av i 

av l _ 

ap 3 

1 ay + *1 

az 

ay 


ax a y 

a *i 



ay 


3*2 3 'l 

_XJT . XZ 


az 


avvg 

'1 ax 



aw 2 

aw i 


aw 3 

+ °2 

-a7 + 0 

3 3x 

+ V, 

ay 



aw_ 



("l W 3 

) + W 2 

IF " 



ap 3 

a 'i 

xz + 

a *2 
Yl 

a *l 

+ 27 

az + 

ax 

by 

3Z 



aw. 


aw. 


+ v. 


+ v„ 


3. Inner layer 
0 ( 1 ) 


0(c) 


U 1 “ V 1 = w, E 0 


p n = P n (x.O.Z) 


au_ 3v„ aw„ 

— ^ + — - + = o 

3x a ; az 
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p 

0(e) 


0 = — 

ay 


* du 2 
f, + — 1 

l _ •» 

*y ay 


0 = — 
3y 


* dw 2 
u + — - 

yz ay 


ax 


o = — 


ay 


fa* 

a^ 3 

ay 

az 


A 

A 

du- 


+ — ^ 

1 2 xy 

3y 


- 0 f A A 

0 = ~ p 3 + { -\ 

ay l yy 


0 = A 

By 


A atV q 

f 2 + ~ 

yz By 


4. Outer- middle matching conditions 

u 7 (x.y.Z) " t/ 1 (x.O.Z) :y - 


30, 


u 2 (x.y,Z) " y -^-(x.O.Z); y 


y 2 


(x.O.Z) + Ug (x , ( 


with similar relationships for w. 


, Z ) : y -» oo 
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V 1 (x.O.Z) = 0 


(x.y.Z) 


au n aw 

~dx 


1 1 Wir n 

I I 1 

- y l-x— + 


az 


(x.O.Z) 


; y - 


£a 2 v, 

v 2 (x.y.Z) ' V 3 (x.O.Z) + ^ (x.O.Z) : y 


(x.Z) = P 1 (x.O.Z) 


3P 

0 = (x.O.Z) 


2 3^P 

P 3 (x.y.Z) ' ^ — j- (x.O.Z) t P 3 (x.y.Z): y - ~ 


t (x.y.Z) ~ 0:y - for n = 1,2,3, •• • 
// 

5. Middle - inner matching conditions 

v 1 (x.O.Z) = v 2 (x.O.Z) = Vg (x.O.Z) = 0 


0 


_a_ 

3y 




A 



(X.y.Z) 


(x.O.Z) 

a 3v 2 

■ y - 3 — (x.O.Z) : y - 
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v 3 (x.y.Z) " y (x.O.Z).y - ~ 

p } (x.Z) = P 1 (x.O.Z) 

P 3 (x.y.Z) " p 3 (x.O.Z) y - a. 

f (x.y.Z) _ f (x.O.Z); y - «. n = 1.2.3. — 
ii n IJ 

u 2 * b 2 ,n y + C 2' U 3 ‘ b 3 ,n y + c 3*‘ K 0 
In y + C 2 , t? 3 " 6 3 In y + C 3 ; y - oo 

with similar relationships for w and w . 

B.C.'s: u n (x.O.Z) = v p (x.O.Z) = t? n (x.O.Z) = 0. n = 2.3. 
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Appendix O: Determination of the 3D logarithmic overlap 

In vector form, we have 


Law -of -the -Wall: — = f (y + ) 

T 


Law -of -the -Wake: \U -U\/u = FCy) 

e j t 


u 


- _ , . dU t df 

From Eqn . ( D 1 ) : — = — 

dY v dy 


. . _ , n „. dU u t dF 

and from Eqn. (0 2): ^ 


Multiplying Eqns. (03) and (D4) by and equating, yields 


t df dF _ x 

Y dy* ‘ ~ V d7 > ~ 


fly ) f (t?) 


where C is a constant vector. 

Rewriting Eqn. (D5) In component form. 


( D1 ) 


( D2) 


(D3) 


( D4) 


( D5) 
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Thus 


rcff dt 

e + — * e 
dy* x dy + 2 


= - 7 ? 


dF v dF 1 

,di) 0 x chj e zJ 


C x e x + C z*z 


L dt 

+ X 

' V 


dF 

V ~ ~ C -* 
dT) x 


f = 


C x in / 


+ 8 


F = 


' C x ,n ^ \ 


( D6) 


dy 


dF z n 
- 71 -d^ = C z 


f z = c z ^ y + B z 


F z ~ ~ c z ln V + 


( D7) 


comparing these results with Eqns. (BIO) and (Bll) yields 
u x = f x = °2 * C x = b 2 ■ B x = 6 2 ■ y+ = y 


(D8) 


F x - - u 2 ‘ A x = ~ c 2' 71 = y 

with similar expressions for the z-components. Note that i] = y dic- 
tates the scale of the outer flow as I = XO. 

Since the Law-of-the-Wall Is 20 (In r^-directlon) . It follows 

that 


+ - 


- 1 . + 

k In y 


+ B 


( D9) 


whence 
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+ . 

COS0 W 

. „ "I 

cos0 w 


= u 

slnK w 

K 

sin*,, 


COS0 , 

sln *w 


B : y - 


comparing this with Eqn. (D8) and its z- counterpart, one gets 


V 


V 


COS^jy 

d 2 


O 

N 


sln *w 


'V 


V 


COS0 W 

* 

l e 2 J 


B z 


s,n *w 


Eqns. (D6> and (07) then Imply that 

-1 » 

F cosy ■ f x " - k cos£ w In 77 + A x :i\ - 0 (DIO) 


-1 » 
- Fsiny s F " - k sintf In 77 + A 

£ nr £ 

If one now assumes that 


1 

A x 

« 

cos^ w 


COS <f> w 

f 

l A zl 

- A 

sln *w 


s,n *w 


then It follows from Eqn. (DIO) that 
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F ~ Yx + " “ K 1 [in V - 2 nj ; 77 - 0 

which establishes, together with (D9). the logarithmic overlap for 
the 3D case. This was also established by the matched asymptotic 
analysis. Eqns. (BIO). (BID and their crosswise counterparts. 
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Appendix E: Details Pertaining to Profile Relations 

El: Determination of the Wake Parameter, n 

Since F 2D Is the defect function of the streamwlse component 
of the skewed boundary layer, the following log. friction law applies 
to It: 


k\ x = In fRe fi /X 1 + kB + 211 


(El-1) 


where 


X x = \/[cos* w ] 1/2 

By using Coles' 111 formulation of the 2D Law-of-the-Wake. subject 
to his normalizing conditions, a further relationship involving n is 

k\ Re _/Re R = 1+11 (El-2) 

x _ o 

°x 

Combining Eqns. (El-1) and (El-2) then yields the equation for n 
In the form 

2n - + K ( s - V] = 0 

6 x 

from which II is determinable, using a root finder algorithm [13]. 


77 


E2: Determination of the function 7(77) 

From the Law-of-the-Wall. Eq. (4). one gets 


dU 


t df 


df v d ,+ 


From the Law-oMhe-Wako. Eq. (7), one gets 


dO _ u r IV dF . rdF _ , dy l 1 

— 1 F 7TZ _x TircosT - F slny-j-M 

dY QU L d7 l } 


Multiplying both expressions by — and equating, yields 


+ df _ 7 ) r fF }dF _ . dyl -1 

y llx - C 0 Sy Jd^ + F siny 


dy 


(E2-1) 


Since the L. H.S. Is a function of y + only and the R. H.S. of 77 only 
(these are independent variables in the context of asymptotic 

-1 


Integration of 


analysis) , both are equal to the same constant, k 
the L. H. S. yields 

f = k 1 lny + <+B 

which is the logarithmic overlap expressed In terms of wall variables. 
Substituting Eq. (6) Into the R. H.S. of Eq. (E2-1) yields 


IL 

U 


2D 


. 2 
Xcos y 


-1 


c "gp F |p laoT dx 

d7 > kcos 2 ? 


= K 


-1 


from which 
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( 9 ) 


dy _ Xcos 2 y \_U [ F 2D _ 1 

d7 > " 4, tanr l K7? Ucos 2 y 


dF 


2D 


c/77 


where U is given by Eq. (7) and F 2 q by Eq. (5). The Initial 
condition required for the solution of Eq. (9) is determined as fol- 
lows: according to the asymptotic analysis the hodograph (Fig. 7) 
starts deviating from the <t> - direction at the point In the logarithmic 

IT 

overlap where the wall- and wake-velocity vectors are In the same 
direction. This point is designated as t\ . 


From Fig. 7, 


+ X 

u = -7 — — 7 

s, n* w 

c ° s *w ♦ 


combining this with Eqs. (6) and (8) yields 


sintf w 

tan0(7?) = 



X 

f 2D (7 » 


(E2-2) 


in the logarithmic overlap. Here 





V 


and i/ + and F 2Q are given by Eqs. (4) and (5) respectively. Eq. 
(E2-2) is solved iteratively for tj q which is the largest 77 for which 
0(ij) = 4> w - y Q is then found from Eq. (8) with V - V Q and 
0 = 4> w - This determines the required Initial conditions for Eq. 
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(9). the solution of which provides y(T}) in the logarithmic overlap 

region, which Is the lower portion of the y(Tj) curve. 

At the upper portion of the defect layer, y = y tEq. (11)]. 

m 

To complete the formulation of y, It Is necessary to devise a 
curve-fit which smoothly connects the lower and upper portions of 
the y(7)) curve. To this end one defines the function 

r s — = i -0 [1-77]“ no) 

m 


The following requirements are necessary to be met by this function: 
1) T(l) = 1 

21 3i>° 


3) llm , = 0 

V - 1 dT) 


4 > ^<0 

dv 


5) " m i ^=0 

71 - 1 d V 2 


(zero curvature as T) -* 1). 

These requirements are met provided that a > 2 and 0 > 0. 
The latter two parameters are determined by matching r to the 
curve described by Eq. (9) at a matching point. 77^. where the 
velocity profile switches from the logarithmic overlap to the wake 
(see Fig. 8) : 
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s v - 


dT_ = -1 dy_ 

dri y mdTj 




where 7(77) is the function at the logarithmic overlap region, and y, 
Is Its value at = 77.,. From this one gets: 


a 


1 - 7 J 


1 cLEL 
1 d1 > 


77=77-, 


and 


p-’.r 

77, is generally a function of Re^ and pressure gradient. Based on 
experimental data, a typical value for It is 0 . 035 . See Appendix E 3 
for an outline of an Iterative procedure to determine 77,. 


E3: Procedure to Determine 77, 

77, Is a function of both Re number and pressure gradient, 
for which a model does not seem to exist at the present time, in 
this work. 77, Is found by an Iterative procedure such that the 
correct value produces a branch of r (77) where r first becomes 
unity at the correct value of 77 for the particular case under Investl- 

I 

gation. In the example of Fig. 11 . the first guess, labelled 77,. 
produces a branch such that r first becomes 1 at an 77 which Is 
lower than the correct value for this particular case (which would 
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typically be taken from an actual experiment). The next guess, 

f I 

V i . leads to a branch In which r first attains unity at an y which 

V t I 

is higher than the true value. Iterating between y^ and 7J 1 even- 
tually produces the correct branch, shown in the figure as a solid 
line. 
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Appendix F: Detailed Analysis of a Juncture Flow 

FI: Preliminary form of juncture equations 

Introducing the normalized variables of Section 2.2.1. (2) 

Into the 3D Navier-Stokes equations, the juncture region equations 
become 


3 U x ref dV X w A x ref 3W 

ax + \ w l ay x^e x^a az 


u 


3U X ref 

ax x w a 


V 


au X w A *ref dlJ 
ay + \ b 6 az 


dP 

ax 


+ 


a 

X w A 

V 

au 

ri/ 

T w 

2 

T 

ax 

*ref 

0 X A 
e w 

ax + 

U 

e J 

XX 


+ 





2 

- 


rU -i 


x ref 3 

v au 

r w 

T 

X w A ay 

U X A 3y 
L e w 

0 

e J 

T yx 


+ 


« *'U, = U/U e . P = P/pO\. T /f = - u] 

x = x/x , . y = y/<x a), z = z"/(x, 8) 
ref w b 


83 



*ref 

a 

[ X * A 

V 

au 

V> 


az 

V 

0 X A 
e w 

az 




\et_ a v av 

V V K V><Jv 4 3Z 

l e w 
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dW *ref w §W X w A X ref w dW _ 

ax + \ w A ay + x & e x^a az 


X w A x ref ap 

V x w A az 


a 

X w A 

v aw , 

rU 

t a 

ay 

x ref 

<5 x a 0X 

0 w 

U 

T w J 



2 

rtf 

T w 

T 


a 

xz 

e J 



*ret 

a 

v aw 

rtf 

T b 

X w A 

ay 

kv ay 

0 

T wl 



2 




e 



t 


X w A *ref a 
x t>° V az 


X w A 


aw 


V 6 o e x w i az 



with the notation 
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V 


€ S 

W 


W 


u 


2 a 

€ € S 

W Hr 


© 


—.2-- C £ ** — 
O.V ' V * V 


© 


U € 

r w " 

these equations become Eqns. (3) of Section 2.2.1. 

Before treating these equations by an asymptotic analysis. It Is 
of Interest to relate C to the parameters e and e. . 

W D 

Both 8 and A Increase with x. At the blade leading edge, 
6/ A " 0. while at the trailing edge 8/ A * 1. For the purpose of 
scaling, the (1/n)th power law growth for turbulent flat plate flow Is 

n-1 

used, according to which 6 ' x n . Thus, for any streamwlse lo- 
cation along the blade. 


* ■ 

This relation may be expressed In terms of X's. using another 
relationship from the (1/n)th power law. namely 

C f - K Re^ 1/n (K- a constant ) 

Thus 
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n-1 

<7 


6 

A 




-E. 


- x 


Re 


x-x 


L. E. 


Re 


x-x 


n-1 

n 


r c x 

' " „ 2( l-n ) 


f b 

- M 

[«*] 

c r 

X. 

e 4 

f 

1 b J 

L w 

w 




and hence 


2( 1-n ) 


. p 

i e » 


l-2n 


2n-l 


Thus £ Is a measure of the streamwise development of the 
corner-flow asymmetry. The case 0/A = 1 formally reduces Eqns. 


(3) to Eqns. (1) and (2), since then also £ s = 1. This 

case corresponds to the symmetric corner flow. In the present 
work, both this and the asymmetrical situtation. which occurs in 
large portions of turbomachinery blade/endwall juncture flows, are 


Investigated. 

The analysis Is carried out for an arbitrary, but non-zero, 
O/A ratio (or C> . so that the leading edge region is excluded from 
the present discussion. 
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The following generalized expansions are defined: 
Middle region 




u ‘ 1 K + e c u 2 + e c U 3 + •*• ] 


V = g [v 1 (x.y.z) + e v + 


c 2 ’ C c'3 


2 

e_v_ + • 


W = f» [w 1 (x.y.z) + e c w 2 + e 2 w 3 + •••] 


P = P 1 (x 'y' z> + V>2 + e c P 3 + 


F // = 


(if 


+ <7 0 . 


f, (x.y.z) + et 0 + e~f. + 


L // 


c 2 // ^ C ° 3 // 


Inner region 


y 

u 

V 

w 


£ F A A 2 A 1 

" f K + e c°2 + e c°3 + * * J 

“ * Pi + € c ; 2 + e c ; 3 + * * * ] 

= f [*, + e c w 2 + £ 2 w 3 + ■•■] 
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A A 2 A 

'■‘"i'V!' Vs + 


T n * 



r € i 




€ I 

. 

L w J 


'iz 


' A A 2 A 

f, + eJ 0 + €_f« + 

l 1 // c2 // c3 // 



Here the parameters f. 


g. h, q. 


t. 


A 

g . 


h. q are In general 


functions of € . e . e . e . C and £ or of some of them; 

C C IT nf 

e and e are small perturbation parameters appropriate to the 
c o 

Juncture, to be determined later from matching conditions; and 


0. 0. ♦ and ♦ are stretching parameters whose purpose Is to 

create symmetric Juncture equations regardless of whether the flow 
Itself Is symmetric or not. These stretching parameters are also 

functions of e„. e„. C and £. The above expansions are 

C C IT IT 

Inserted into Eqns. (3) . yielding: 

Middle region 


f -f- (u n +♦ • *) + 
3x 1 


C SL (y J. » . . ) 4. — — ( W t • 

5 0 3y V 1 + ’ + c 9 az 1 


•) = 0 


(U, +••♦) <u 1 + | (v 1 


1 


) — (u. + •••) + 
6y 1 


HA ( W + 
c* 1 


> tz (U 1 + 


) = 


- £ <p i + 


> + 


a [f 2 a a ... 

3x [£ e w € w 3x <u 1 


O + e c f, + *• 

1 xx 


89 


f a ff 2* a , 2 r 

fa? I? e w e wTy (u l + - ) + e c f l + ••• 

1 l yx 


. _i_ a f 2 a a , 

+ c* az c* 6 w e * dz (u i 


+ + C ®( fl zx + 


f0(o i + ...) ± (Vl t •••> + « ^ (v, + ...)| 


Ml +■ A < Vl + ...) = -f ± ( Pl 


8 0 2 A d p r 

+ dx If e w € * aj (V 1 + "* ) + e c f i + * “ 
1 l *y 

£_ 3 fg 2 * 3 , 2 r 

* ay [* € w e w ay <v i + **° + € c f f i + *** 

L l yy 

{ 0 f fl 2 A 3 Pf 

* « ai [c* E * E ». ai <v i + ' ’ ’> + E c f'l ♦ 

f/lCO, + •••) ^ Of, 


-£*- "l £ <*! » •••> - - ££ <P, 


♦4 


/> 2 ^ a, 

? e * 6 * a? (Mr i + 


•) ♦ C 2a e; 


M- 


H- 

"H 

; (y 1 + • • •) + 
• • • ) + 

ll* 

]]• 

•••]] 
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ii i e 2a ± (w 
<f> ay <t> w*w 3y 1 


+ 


, c- 4 




2 * a , 
£ *<V a; ‘"l 




imposing symmetry In Y and Z and In V and W on these 
equations results In the following relations: 

2 



f = -\. g = h = e c 

* = - « e c 


and 
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f = 


Inner region 


2- 

e w € >y 

2a 

£ c € c 


f dx ( °1 + "* ) +« * “ + -**>+ 1 7 ji (*,+••• 

ay zt oz 1 


^ «j A A 

f < "l + *** ) If <"i + •••)+€ fif +...) A (u 

0 3y 


1*1 

C* 


~ ( W + •••) -2; (O +...)=- f- (p, + 

♦ 1 3z 1 3x1 


. a i 2 a a a 
3x £ € w € iv 3x + 


V 2 a 

■> ’ 4 'i 


K, * •••)] 


1 

$ 


^ K ^ ty 

( f V £ ■■■]] £ 


+ • • ♦) + 


2 

s c n 


£ 3 f 1 2 a a A 2 a 

„JA A A e * € * ~ a <w i + *•*> + (f + 

C* 3z Ic^ 3z 1 c 


aaa n a2 

f° (o 1+ ...) A CV, +...,+« (v 1+ ...) (V + 

? 3y 


) =0 


+ •••> + 


•) + 



• • • ) + 
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(w, + •••) — (v, + *••) = 


a 7 

C* 


V » T 

a"' 1 

3 z 


£ 3 .^ 

- jL — (p + • * • ) + 

* * i 

0 ay 


A fi e 2 * A ( ; 

ax [* e w e w ax i 


■> + 4 h + 

1 xy 


IV "H 


L A 

A A 

0 ay 


g 2 * 3 . A 

a e e — (v, + 
* w w * 7 

0 ay 


) + 


£ * lv -)] 


£ a a 2 * 3 / A 

1 — * -®- € € ( V- +•••)+ 

A A A W W A 1 

t* az U* * w 3 z 


IV-|] 


^ . * . a / * 

fh (U, + • • ♦) -r- (W. + 
1 3 x 1 


■) + £ (v^ + •••) •— ( +•••) + 


0 


ay 


aP 

*/» * a , * f a , a 

•* — (w. + • * • ) — (w. +■•-)=- -r— (p. +•••)+ 

A 7 A I A 02 7 

C* az c* 


ax 


h 2* a . * . t 2o 2 

TV»5 i, i , '"" ! £ w 


‘ 

2 

r e c 1 

A 

— 

+ Q 

KJ 


*i * • 

xz 


II 
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i. 

$ 


d f h 2 * 6 * 

h e w € w (w i 

0y Qy 


> ♦ < 2CT 4- 


2 I 

f e c] 

r* 1 

— + q 

f + ••• 

U e wJ j 

l n ] 


iili € 2 a A 
ci a$U* w w 32 1 


+ — ) + C 2a e 2 - 


Si H['V-] 


Imposing symmetry again, one obtains 

A A „ A A AAA A 

<7=3. f - 1. g=ft = c c 6 c . 0 = C* « €c c e c 

With the scales thus determined, the equations are rewritten 
as follows: 

Middle region 

ai ( U 1 + e c U 2 + * * * ] + Ji ( V 1 + e c V 2 + * * * ) * 

Tz h + c c *2 + * • * ) = 0 



3 r 1 

(“l + € c °2 + •••] 

a; (“l + e c “2 * • • • ) 
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te c *2 ■'•••) 4} ("l +e c u 2 t "-) + h ,€ c "V']' 
("l + £ o “2 + • • • ] ■ ■ S ( P 1 + € o p 2 + • • • ] f 

Tx [4h i("l + Vs + ••• ) * £ c('l„ + € o' 2 „ ’ ] l 

'iMh + e c “2 + •••) ♦ e o (%, + *c\ * ]] ♦ 


fz a c ai h + Vs + •••) + «. \\ + Vs + 


c c[["i t£ c u 2 + --) ix h +t o'* + "-) t h t v* + -"] 


a? <V 1 * Vs + •••’ + 


[*, <! [ * ! *-)T:(VV! , ''']l = 

- -Ty [' , l*V! , -") , l;[Vcl('! ,i o's t -) t 


lv ! -]l 


a 2 a a , 

+ — € € -r - 1 (V n + 

6y c c 3y 1 


0 + «c \ ♦ 
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O fo 


il 


ty <"l + 


* • ) + 


0 f 2 a 9 t Or 

3z[ e o £ o Tz * 

[[“l + •••] T (w 1 + •••> + <»,+ 

<"l + •••> Tz <w l + •••>] * 

"^[ f ’i 1 ' e c p 2 t6 o p 3 + '") f |;[ £ c ? c ~h ‘"’i + •••> + 

*«h« * ■••]! ♦ 

+ ■§?{'% Ty lw 1 * •••> * 4 |', y2 * ••■]] ♦ 


, a 2 a 3 , 

+ -TZ e e — (w . 
az c C dz 1 


Inner region 


* ---» * * ---]] 


dx K + € C S 2 + • • * ] t ^ (v, + e c v 2 + • • • J + 

^(*1 + e c ^ 2 + •••]= 0 

02 


A A 

e c «/, t • 


\ 0 , A 
O S (U, + 


•) + (V 1 + 


•> ^ c, + 

3 y 


) + 


(^ 1 + 


, a a 

7 - i * 

02 


]■ 
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* 6 , A . 3 I 2*2 6 , * 

" *cTx (p i + '" )+ a^l € c € c a? (t/ i + ‘** ) + 


2* f* 

€ o € c 'l ♦ 

. XX 


0 0 a 

tE c 


ft, * ...1 

[ yx ) i 


a*U* 1 1 


♦ c o"a T e c“3 f • 


f MV 


V2 » f 


• ) + + 


0 * 

•) ( V- + • • • ) + 

* 1 
ay 


+ 


) — (V, + 

a A 1 

8z 


A 2 * 

V>2 + € c P 3 


e ' e ‘ I'v ’ 


-]] 


0 I 2* 0 a 2 f * 

74 e c e cT A % + •••> + e c h + 

ay 1 ay l yy 
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±\*% ± <v\ 
8zl ° C 8^ 1 


) + €‘ 




Iv.l’l 


(u. + •••) ■§- (W, +•••)+( V. + • • ♦) — (*,+•••) + 

1 ax 1 1 8y 1 


, A 3 , A 

( w . + • • •) — <w, + • • * ) 
1 « 

3z 


3 f A A 2 A 

— l p l + e c p 2 + e c p 3 + 


0Z 




4 a 3 8 . A 

e„e — (w. + 
c c 8x 1 


) + 


'"U 


8 2 - 8 . - 
— G „ G — (w., + * 
-I c c - 1 

8y l 8y 


8 f 2- 3 . A 

— 1 e € — (wr, 

8zl c c a? 1 


) + e 


2 f A 


-]] 


8zi dz 

Sorting the above equations according to orders of magnitude 
results in the following equations: 

Middle region 
0(1) 


8(/ 1 8v. 3tv 

■aT + -aF + -aT = 0 
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3u l 

aw i 

3u l „ 

3p 1 

17 + *i 

dy + W 1 

dz 

" 17 


du 2 dv 2 dw 2 
dx + dy * dz 


a 2 

— (u,u 0 ) + V- — — + v 0 
3x 1 2 1 dy 2 


0f i a 'i 

_Je _Jk 

dx dy dz 


at/ 3 ay 3 awr 3 
dx + dy + dz 


Tx (0 l t/ 3 ) + U 2 


du 2 dUg 

— r — + v. — — + v_ 
dx 1 8y 2 
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3u 3 at/p 3t/ 1 

w -r^ + w 0 — ^ + W 0 -r— 

1 8z 2 3z 3 3z 


ff3 

ax 


at. 


3f, 


at. 


XX 


ax 


u i ax + v i 


a *i 

w i “ST + v i 


ay 

3v 




zx 


az 


av. 


ap. 


at. 


at. 


IT + "i 

3z 

ay + 

ay + “ 

3z 

a "i + 

-5T + w i 


ap 3 

a 'i 

Y z * 

a 'i 

zz 

3z 

az + 

ay 

3 z 


Inner region 
0(1) 


A 

d</ 1 

A 

fri t 

A 

3 "1 

ax + 

A 

ay 

az 


= 0 


2 a 2 a 

an. a^u, 

o = t 


- *2 
3y 


dz 


(FI) 


0 = 


ffl 

3y 


0 = 


3p 1 

3z 


0(e) 

c 
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CVJ o 


3o_ 3Vp dw p 

* t — - t — - = 0 


3x 


ay 


dz 


0~3- 

ay 


du 2 * 

— £+r n 
* 1 
3y yx 


az 


A 

di/ 


2 * 

tr. 


az 


zx 


0 = 


0 = 


apg 

ay 

A 

apg 

9z 


9w_ 0v_ aw* 

+ — + = 0 

djr a a 

9y 3z 


0 =— 

ay 


a “3 * 

a 

a “3 f 

— +'j> 

+— 

~ +t 2 

l ay yx J 

dz 

l az zxj 


ap 8r i af i 
JCL + 2 


o + - 

A A 

3y 3y 


dz 


a P * 3, i a 'i 

_yz zz 


0 = + 

A A 

dz dy 


dz 


F2: Matching procedure 
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The following relationships exist between the juncture coordi- 


nates and those of the endwall and blade: 


„ 2o_ y Vi, 0 L y _t 


Z c ^ _ c Z~ £ X b 0 g b € b z b 
c 'V‘ f V £ c' { X b 6g c 


A 

« e c b 


X w A€ w € w y w g w g w 


y c - ^A«e c 


{£ c Y W 


iVf, 

- { k^Oc, - £c 0 


The first pair of relations Indicates that for smooth matching 
between the Juncture and the blade. 


Ilm y = 


llm z = 


c c -* e b (blade matching ) (F2-A) 


The second pair indicates the requirement for smooth match- 


ing between juncture and endwall: 
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lim y = 
y -*0 


llm z = 


e wV 

Z 


w 


'w 


e -» —r ( endwall matching ) (F2-B) 

C 6 


A re/ 

Since | s - — r- and from Eqns. (F2) It follows that 
*w 


llm { = !-* llm e c = e w 

Z -*00 Z -«30 

c c 


lim 4 = C -* lim e = e. 


(F3) 


In a similar manner 


llm 


A 

z 


€ c = V 


Eqns. (F3) are 


Mm e c 


A 



y c ^o 

in agreement with the previously obtained 


2 a 

e e 

result 4 = 

€ c € c 

Thus € and e have the correct limiting values to guarantee 
c c 

smooth matching with Regions I on both endwall and blade. Hence, 


one may define e and e In a similar manner to the corresponding 
c c 

A A 

definitions of e^. e^, e b and e^. namely 
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u 

r 


c A V 



where u and A are. respectively, friction velocity and length 
T c c 

scale characteristic of the juncture region. 

Since the equations of motion derived earlier are symmetric. It 
suffices to perform the matching procedure for one half of the Junc- 
ture region. The matching will be done with the endwalf. along the 
following limiting lines: 

1. Juncture middle region and endwall middle layer; 

2. Juncture inner region and endwall Inner layer; 

3. Juncture middle region and endwall Inner layer. 

These lines are denoted by the numerals 1. 2 and 3 in Fig. 

4. A completely analogous procedure applies for the blade- 
matching, denoted by the numerals 4. 5 and 6 In the figure. 

In the following, the subscript w is dropped from e^, since 
both juncture- and endwall- expansions are using this same pertur- 
bation parameter, as demonstrated earlier CEqns. (F4)j. The ex- 
pansions for the endwall layers Incorporate the results found earlier 

for the case of pressure-driven boundary layer (see Appendix C) . 

Y 


2 
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1 . 


Z -* 0 
w 


Z c 


U w = U W <x -y- z w ) + €u w i ‘ eZu w +•** = ^ <x.y.O) + 

1 w w 2 *3 w -\ 


3w, 


wr 


2 

2 a u 


1 2 Z c W 1 

-37— (x,y.O) + € ^ 5- <x.y.O) + 

W 


az 


w 


• + €u <x , y . 0) + 

W 2 


du w 

2 2 2 
« z e -37“ <x.y, 0 ) + --- + e w^ (x.y.O) + 

tv 3 


w Ci (x.y.z c ) + 


€U (x.y.z ) + e u (x.y.z ) + — ; z ~ 


u (.x.y.z ) 
C 1 


U w 1 <x • y ■ 0) ; 2 C - 00 


a«/ 


w , 


V ‘ 2 c l2~ *c - 


mt 


% "'-V 


2 

2 a u Bu 

Z C *1 W 2 

-5- «“ (x.y.O ) \z (x.y.O) + 

c dZ tv 


w (x.y.O); z - 
w 3 0 
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(x , y , 0 ) + 


V = ev t e v + 
w w 1 w 2 


av » 

= « <*./.<» + «% -jjj 

I IT 


+ € v +••• - ev^ (x.y.z c ) + e v (x.y.z c > + 


: z 


v (x ,y . z ) 
C 1 


v ix.y.Q): z -* «» 
7 c 


0V 


V C 2 <*' l '' z c > 


az^ ‘ X - I '- 0> *V ! : Z C 


W w - "V, ”v, te ’V 3 + 


ew „ +€ w + 
C 1 C 2 


2 -*0 
w 


Z -K» 

c 


Differentiate w. r. t. Z : 

w 


aw . 


W , 


a 2 w 2 2 a 3 w 

1 "i e Z C w i 

-rz — (x.y.O) +ez c (x.y.O) +— 7 ; g— (x.y.O) + • 

dZ^ 0Z 

VV IV 


04 V 

W 2 2 
+ € -7= — (x.y.0) + c z 
az * c 


d 2 w 


w . 


(x.y.O) + * 


3 Z 


wr 


p ”o 

• + e -r ^- 2 (x.y.O) + 
3Z w 


a*v aw c 3w c 

al “ tx.y.z^jJ + c-^ (x.y.z c ) + € 2 -g~ (x.y.O) + •• -z c - - 
c c c 
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dw_ 


aw. 


w-, 




-jj- <x.y.o>: x. -- 


dw 


d 2 W 


W n 


0W 




-r— ^ ( x.y.z ) 
3z_ 7 c 


3Z 


(x,y,0)+— — (x,y,0); z c - 


w 


w 


aw 




Q 9 

2 a w w aw w 

■ 4 5T (x.y. 0 )+ 2 r ^(x.y. 0) + 

^ az c az 

w w 


dw. 


w„ 


-=■= — <x , y . 0 ) ; •* co 

dZ w c 


3P 1 

P = p i (x.0.0)+ez^ -r-r- 1 - (x ,0,0) + 
w i c az^ 


2 2 

2 Z c 3 P 1 , o o> 
e -g 2 ~ ( x,0,0) + • 


az 


w 


2 

t e p (x,y,0> + *** ~ p (x)+ep (x) + 
w 3 C 1 c 2 

€ 2 p c (x.y.z c ) t z c - ~ 

3 

p (x.y.z ) - P. (x.0.0) 

C. C I 


(F4) 
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(F 5 ) 


3P 1 

p c 2 <x ' y ' z o > " z c az^ <x '°' 0); z c - 


Pc 3 (x.y.z c ) 


Z 2 a 2 p 

~ (x. 0 . 0 )+p ty (x.y.O): z c 

dZ w 3 


W 


■t (x.y.O) + cz c — — 

it w li c dZ w 


a V, 2 a V 

1 (x,y.O> + € ^ (x ,y ,0) + 


2 2 
3Z 


a s 

+ €f w (x.y.O) + € 2 z c — 07 ^ (x.y.O) + 

2 // * 


e (x.y.O) + 

% 


f (x . y . z„ ) + et (x . y . z ) + 
c. c c_ c 

// // 


€ f c q <xy 'V + 

3 // 


: z c - « 


f Cl (x ' y 'V V (x ' y0); 2 c ** 

// 7/ 


3f. 


w. 


f (x.y.z c ) 
// 


'// 


2 c"lz"" (x.y.O) +f w (x.y.O): z c - 




7/ 


f c_ (xyz c ) 

3 il 


- _C V/ 

2 2 

az 

w 


at. 




(x.y.O) +z 


7/ 


■e az <*-y.°>+*_ '*■> 

* 3 (/ 


, 0 ) ; z •* oo 
c 
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2 -*oo 

C 


U = eu (x.y.O) +0(e) +e 2 t/ (x.y.O) + 

w w_ W 3 


A > A A x i 

'"o, <x > ,z c > + 


t (X./.V + -- 


* 7 -4 00 

• c 


u (x.y.z ) "0; z -* « 
c c 


£ (x.y.z ) - (X.y.O): z 


A A 


c W 


£ (x.y.O " w (x.y.O): z 


A A 


c w 


2. . _3' 


= e ev (x.y.O) +0(e ) te ev (x.y.O)f- 

U/ W- 


ee£ (x.y.z) + e 2 ev (x.y.z ) + eev (x.y,z c ) + 


Z -* oo 
C 
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( F7) 


A , A A A 

V., <X.)-. z 0 ) o: z c - 


A A A 

% ">'V 


\ '* v.O): z c 


A A A 

v (x . y . z ) 
Co c 


V, <x >’- 0): 2 c 


w * ■ £ “V <*•»•*„> + ** *_ 

^ 3 


e ®^c, <x 'y^ c ) + e 2 e w c (x.y.z c ) +e 3 evP c 

2 3 


'1 


Z -*0 
w 


Differentiate w. r. t. Z 


tv 






2 *3 


— <x .y . 0) +0 <e) + c —Q 2 — (x.y.O)- 

W IV 


aw„ aw 

C, C- 

1 , A A 2 . A A 

— — ^.y.z c ) + e * <*-y-v + 

3z c az c 


no 


- A A 

<x , y . z ) + 
c 



oo 


a* c 

2 3 , A A A 

6 __ <x.y.z c >: z c - 

3z c 


aw 


a z 


1 A A 

- <x,y,z ) 
c 


0 ; 


dw„ 


aw.. 


w. 


a 


(x . y . z c ) az 


(x.y.O) : z - °° 


w 


aw 


aw, 


az 


3 . A A 

— (x.y.z ) 

A c 


w. 


az 


(x.y.O): z„ - 


w 


P - P y (x , 0 . 0 ) + O ( e) t e^p (x.y . 0 ) + 
w i w 3 


p (X ) t 
C 1 


ep (x) + 
C 2 


P A A A 

+ c p c (x.y.z > + 
c 3 ° 


• : z_ — 


p (x.y.z ) = P, (x. 0 . 0 ) 
c . c * 


( F8) 


( F 9 ) 
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A A A 


A A 


P w (x . y . 0) : z 
*3 c 


r w = (x.y.O) +0(e) +€f (x.y.O) + 

ii 1.. "p 


U 


2 a a 

e *, (x.y.O) + 

% 


, «. f . A A t A AA 

+ *** t (x.y.z )+ef (x.y.z ) + 
C 1 c c o c 

// 2 // 


2 A a a 

e * Cx.y.2 ) + •..; z - o» 
c c 

3 // 


A y A A * «. A A A 

f, (x.y.z c ) f w (X.y.O): z„ - 




'll 


A . A A s - A A A 

**„ t w (x.y.O): z„ - 


"// 


V/ 


A A A 


t n (x.y.z c ) - r w (x.y.O); z^ - 


// 


3 


H 


(F10) 
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y c -° 


°c <x ' y ' z c ) + eu c <x,y,z c ) + € u c (x *y* z c ) + 
12 3 


A A 2 A A I 7 ” 

' e u (x.y.Z >+e u (x.y.Z ) + •••: j 

w 2 *V Wg W ly -*co 


Differentiate w. r. t. y: 


)u 3u„ 3u 

®1 c o 2 C 3 

— — (x.y.z ) + e — r — (x.y.z ) + e — — (x.y.z ) + 
3y c 3y c 3y c 


a v a v 

■f — — <x.y.Z w )t€ — — (x ,y .2^) + 

6 L 3y 3y 


du w 

2 4 a 

€ — r <x.y.z w > + 
3y 


In analogy with the* discussion in Appendix B. 


"ay (x 0 z c ) " 0 


9 ° c 2 . at/ °3 . V x - Z c ) 

3y y ' ay y 


; y -♦ o 


113 





ay 


ay 


• • • * y -* 


or 


« c ’ 0 / (x./ c ) In y-l-^ (x.z c ). / = 2.3. •••; y-0 
% ’ */ (x - Z w ) ln *+*/ (x Z w ) - 1 = 2.3. •••; y-« 


V c = (x.y.z c )+e <x.y.z c > + 


2 A A A 

e ev (x.y.Z 
*2 " 


3a a 

t € € V 


» 9 + { 


y-*0 

A 

y -*oo 


Differentiate w. r. t. y: 


av 3v 

C 1 2 °2 
e -3- <x.y.z c ) + e —u.y.z c ) + 


dv w 

2 2 * 

* e (x.y.Z ) 

- w 


By 


t e 


Bv 

3 w 3 


By 


< x -y - z w > + 


y-0 

A 

ly -k» 


3 v 


IT (x 0 z c ) = 0 


3v c 


Bv. 


Wr 


ay 


(x.y.Z^) : y - - 


(Fll) 


) + 
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W c as e w Q (x , y . z c ) + e w c <x . y . z Q ) + • • 


" e (x.y.Z^> + 


+ € 


{ y-*0 
* 

y-«o 


w (x.O.z ) ' w (x.y.Z >; y - 
C/ c w /+1 " 


~. I = 1.2. 


P = P, (x . 0. 0) + ep (x) + e p (x.y.z ) + 
cl c 2 C 3 


ry-»0 


P^x.O.Z^ + c p w 3 + 1;. 


p =0 
C 2 


( F12> 


P C 0 (X '°' 2 C ) ' P W , <Xy 'V ; y - 

3 o 


r =r (x.y.z > + er (x.y.z ) + ef (x.y.z ) + 


°ll C 1 


// 


"// 


>11 


L <*•#•*„> +« 'V + f 

’// // « 
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"// 


(x.O.z ) 
c 


A 

'w 


(x.y.Z^) ; y - 


n = 1,2. 


// 


Important results from matching procedure: 

1. From Eqns. (FI). (F6) and the one equivalent to 

(F6) for y - oe> plus the boundary conditions u = 08 

c c. 


z =0 
c 


or 


y c =0 


It follows that u Is described by Laplace's equation with 
C 1 

homogeneous boundary conditions. By the Maximum and 

Minimum Theoren. u m o. 

C 1 

2. From Eqns. (F4) and (F9). both p and p are 

C 1 C 1 

given by the external pressure evaluated at the corner. 

Moreover, since both p and p vanish [ Eqns. <F10) 

°2 C 2 

and (F12)]. the above conclusion carries through order e. 

3. By Eqns. (F5) and (F12). dP^/dZ w vanishes along 
the corner. From symmetry. 0P^/ay^ also vanishes along 
the corner. 

4. From Eqns. (F11) and their equivalents on the blade, 
the streamwlse velocity has a common logarithmic overlap 
between the corner middle region and the Inner layers on 
both endwall and blade. 
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Appendix G. Effects of Rotation 


«** 

Y 



If the (x.Y.Z) coordinate system Is rigidly attached to a ro- 
tor. rotating at a constant angular velocity fl around the x-axis. the 
inertia terms of the N-S equation will include centrifugal and Coriolis 
accelerations and will assume the following form [14]: 


O^ + y^ +l y^. n 2 fl -0fl 


6x 


ay 


az 


ax 


u — + V— +IV— -2nw -n 2 fl- 


ax 


ay 


az 


ay 
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0 m t v * an V - (? R S& 

ax ay az' az' 

After the usual normalization, these become 


(J 2 f ad ] 2 

^ [u su . v m t w rn . 1 2faJ_ „ is 

XA l 3x 3 y 3Z J XA 3x 


(J 2 


t-t ft/ “ + V 
XA l 3x 


gf + «-Sf)-.nd. « 


Ki!„ an 
xa ay 


<J 2 


-fi fuM +v aw + w Ml + 2nu - v _ 
xa l ax ay az J + 44 11 u e 


K ) 2 „ an 


XA 


az 


_p 

Multiplying through by \6/U , as done before for the non- 

o 

rotating case, yields: 


u it/ au au 

ax ay az 


n ft. 


u 


aw 

ax 


v |^ +w |^_ 2 nxA w _ 

ox dY dZ y 

e 


nft. 


17 


aft 

ay 
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2 


u m + v w m + k v 
ax ay az ^ 

6 


In turbomachinery components. 



R 


M 

az 


OR 


0 


0 

e 


0 ( 1 ) 


Also. R »A and hence R scales with the outer, Invlscid scale: 
o o 


R = O(XA) - — = 0(1) 

° 0 

e 

With these non-dimensional parameters thus established, it Is 
now possible to determine the effect of rotation on the various sub- 
layers and sub-regions found earlier. This is done by the usual 


2 

sorting out of terms In each sub-region to orders 1 . € and e . 

The following momentum equations are thus found for the case 
of a rotor, rotating at constant angular velocity about the x- 
coordlnate: 

1 . Outer flow 


U 

U 


Bx dY BZ dx dx 


17 * v §7 * w § - 2 w - n I? 


ap 

ay 
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2. Middle layer (Region I): 
0 ( 1 ) 


- ft 


aft 

az 


3ft 

az 


0U 1 . „ 

a<i i 

8u i 


a/? 

ap 1 

0X 1 

17 + "i 

az 

- rt 

0X 

“ 17 




o 

dR 

ffi 




n 

dy " 

ay 


a "i + 

dw -\ 

- rt 

aw 

3p 1 

“T + V- 

0X 1 

17 + "i 

az 

0Z 

” IF 


0 ( 6 ) 


o(e 2 ) 


Tx <U lV + v i 


du 2 

01/ ^ 

3w 2 

du l 

17 + V 2 

-r — + W. 

0y 1 

IF + w 2 

az 


apg 

dx 






ay 


ap 2 

2 *1 = 17 


aw_ 

U, ~Z~ + Pr 
1 8x 2 


dw ^ 

a^ 2 

a *i 

”~r — t V, 
0X 1 

17 f v 2 

ay 
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3. Inner layer (Region I): 
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4. Middle region (Region II) 
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5. Inner region (Region II) 
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From these equations, 
served: 


at. 


22 


A 

dz 

the following 


effects of rotation are ob- 


1. In the outer flow, there Is a centrifugal term In the 
streamwlse component of momentum and both Coriolis and 
centrifugal terms In the crosswise and normal components. 

2. In the middle layer of Region I. there are centrifugal 
terms to 0(1) In all three momentum components, and 
Coriolis terms In the normal- and crosswise- components to 


p 

orders e and e . The streamwlse momentum component Is 
not affected by rotation to these orders. 

3. The Inner layer of Region 1 Is unaffected by rotation 
except for the 0(1) normal momentum, where the pressure 
gradient Is balanced by a centrifugal acceleration. 

4. In the middle region of Region II. there are centrifugal 
terms in all three momentum components to 0(1). In 
particular, crosswise pressure gradients are balanced by 
centrifugal accelerations. There are no rotational effects 
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on the 0(e) equations, while Coriolis terms reappear In 

p 

the two crosswise momentum components of O ( e ) . 

5. The only effect of rotation on the Inner region of Re- 
gion II Is In the 0(1) crosswise momentum components, 
where pressure gradients and centrifugal accelerations bal- 
ance each other. 


126 




1. Report No. 2. Government Accession No. 

NASA CR-3792 

3. Recipient's Catalog Nc. 

4. Title and Subtitle 

Scaling and Modeling of Three-Dimensional, End-Wall, 
Turbulent Boundary Layers 

5. Report Date 

June 1984 

6. Performing Organization Code 

7. Author(s) 

Uriel C. Goldberg and Eli Reshotko 

8. Performing Organization Report No. 

None 

10. Work Unit No. 

9. Performing Organization Name and Address 

Case Western Reserve University 

Dept, of Mechanical and Aerospace Engineering 

Cleveland, Ohio 44106 

11. Contract or Grant No. 

NAG3-270 

13. Type of Report and Period Covered 

Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

505-31-0A (E-2024) 

15. Supplementary Notes 

Final report. Project Manager, Peter M. Sockol , Fluid Mechanics and Instrumenta- 
tion Division, NASA Lewis Research Center, Cleveland, Ohio 44135. This report 
was submitted by Uriel C. Goldberg as a dissertation in partial fulfillment of the 
requirements for the degree Doctor of Philosophy to Case Western Reserve Univer- 
sity, Cleveland, Ohio. 

16. Abstract 

The method of matched asymptotic expansions is employed to identify the various 
sub-regions in three-dimensional, turbomachinery end-wall turbulent boundary 
layers, and to determine the proper scaling of these regions. The two parts of 
the b.l. investigated are the 3D pressure-driven part over the endwall, and the 
3D part located at the blade/ end-wall juncture. Based on the results of the 
analysis for the pressure-driven part, models are proposed for the 3D Law-of-the- 
Wall and Law-of-the-Wake. Comparisons between these models and the data of 
van den Berg and Elsenaar and of Mueller show good agreement between models and 
experiments. 

17. Key Word* (Suggested by Author**)) 

Boundary layers 

Turbulence 

Turbomachinery 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 34 

10. Security Classlf. (of thl* report) 20. Security CImsII. (of this page) 21. No. of pages 22. Price* 

Unclassified Unclassified 134 A07 


For sale by the National Technical Information Service, Springfield, Virginia 22161 

NASA-Langley , 1984 






National Aeronautics and 
Space Administration 


THIRD-CLASS BULK RATE 


Washington, D.C. 
20546 

Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




